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We describe a recognition algorithm for a subset of binary linear context-free rewriting systems 
(LCFRS) with running time where M{m) = 0{m^) is the running time for mx m 

matrix multiplication and d is the "contact rank" of the LCFRS - the maximal number of 
combination and non-combination points that appear in the grammar rules. We also show that 
this algorithm can be used as a subroutine to get a recognition algorithm for general binary 
LCFRS with running time The currently best known oj is smaller than 2.38. Our 

result provides another proof for the best known result for parsing mildly context sensitive 
formalisms such as combinatory categorial grammars, head grammars, linear indexed grammars, 
and tree-adjoining grammars, which can be parsed in time It also shows that inversion 

transduction grammars can be parsed in time In addition, binary LCFRS subsumes 

many other formalisms and types of grammars, for some of which we also improve the asymptotic 
complexity of parsing. 

1 . Introduction 

The problem of grammar recognition is a decision problem of determining whether 
a string belongs to a language induced by a grammar. For context-free grammars, 
recognition can be done using parsing algorithms such as the CKY algorithm (Kasami 
1965; Younger 1967; Cocke and Schwartz 1970) or the Earley algorithm (Earley 1970). 
The asymptotic complexity of these chart parsing algorithms is cubic in the length of 
the sentence. 

In a major breakthrough. Valiant (1975) showed that context-free grammar recogni- 
fion is no more complex than Boolean matrix multiplication for a matrix of size m x m 
where m is linear in fhe length of fhe senfence, n. With current state-of-the-art results in 
matrix multiplication, this means that CFG recognition can be done with an asymptotic 
complexity of 0(n^-^®). 

In fhis paper, we show fhat fhe problem of linear confexf-free rewriting system 
recognition can also be reduced to Boolean matrix multiplication. Current chart parsing 
algorithms for binary LCFRS have an asympfofic complexify of where / is fhe 

maximal fan-out of fhe grammar.^ Our algorithm takes time for a consfant d 
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1 Without placing a bound on /, the problem of recognition of LCERS languages is NP-hard (Satta 1992). 
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which is a function of the grammar (and not the input string), and where the complexity 
of n X n matrix multiplication is M (n) = 0{n‘^). The parameter d can be as small as /, 
meaning that we reduce parsing complexity from 0{n^^ ) to 0{n‘^^ ), and that, in general, 
the savings in the exponent is larger for more complex grammars. 

LCFRS is a broad family of grammars. As such, we are able to support the findings 
of Rajasekaran and Yooseph (1998), who showed that tree-adjoining grammar recogni¬ 
tion can be done in time 0{M (n^)) = 0(n‘*’^®) (TAG can be reduced to LCFRS with d = 
2). As a result, combinatory categorial grammars, head grammars and linear indexed 
grammars can be recognized in time 0{M{n?)). In addition, we show that inversion 
transduction grammars (Wu 1997, ITGs) can be parsed in time 0(nM{v?)) = 0(n®'^®), 
improving the best asymptotic complexity previously known for ITGs. 

Matrix Multiplication State of the Art. Our algorithm reduces the problem of LCFRS 
parsing to Boolean matrix multiplication. Let M (n) be the complexity of multiplying 
two such n X n matrices. These matrices can be naively multiplied in O(n^) time by 
computing for each output cell the dot product between the corresponding row and 
column in the input matrices (each such product is an 0(n) operation). Strassen (1969) 
discovered a way to do the same multiplication in 0(n^ ®^°^) time - his algorithm is a 
divide and conquer algorithm that eventually uses only 7 operations (instead of 8) to 
multiply 2x2 matrices. 

With this discovery, there have been many attempts to further reduce the com¬ 
plexity of matrix multiplication, relying on principles similar to Strassen's method: a 
reduction in the number of operations it takes to multiply sub-matrices of the original 
matrices to be multiplied. Coppersmith and Winograd (1987) discovered an algorithm 
that has the asymptotic complexity of Others have slightly improved their 

algorithm, and currently there is an algorithm for matrix multiplication with M (n) = 
0{n‘^) such that a; = 2.3728639 (Le Gall 2014). It is known that M(n) = fl{rf logn) (Raz 
2002 ). 

While the asymptotically best matrix multiplication algorithms have large constant 
factors lurking in the O-notation, Strassen's algorithm does not, and is widely used in 
practice. Benedi and Sanchez (2007) show speed improvement when parsing natural 
language sentences using Strassen's algorithm as the matrix multiplication subroutine 
for Valiant's algorithm for CFG parsing. This indicates that similar speed-ups may be 
possible in practice using our algorithm for LCFRS parsing. 

Main Result. Our main result is a matrix multiplication algorithm for unbalanced, 
single-initial binary LCFRS with asymptotic complexity M{n‘^) = where d is 

the maximal number of combination points in all grammar rules. The constant d can be 
easily determined from the grammar at hand: 

( ip{A)+ip{B) -(p(C), 'I 
d= max max < ip{A) — (p{B) + ip{C), >. 

{-^{A) + ^{B) + ^{C)\ 

where A ^ B C ranges over rules in the grammar and ip{A) is the fan-out of nonter¬ 
minal A. Single-initial grammars are defined in §2, and include common formalisms 
such as tree-adjoining grammars. Any LCFRS can be converted to single-initial form by 
increasing its fan-out by at most one. The notion of unbalanced grammars is introduced 
in §4.4, and it is a condition on the set of LCFRS grammar rules that is satisfied with 
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many practical grammars. In cases where the grammar is balanced, our algorithm can 
be used as a sub-routine so that it parses the binary LCFRS in time A similar 

procedure was applied by Nakanishi et al. (1998) for multiple component context-free 
grammars. See more discussion of this in §7.5. 

Our results focus on the asymptotic complexity as a function of string length. We do 
not give explicit grammar constants. For other work that focuses on reducing the gram¬ 
mar constant in parsing, see for example Eisner and Satfa (1999), Dunlop, Bodenstab, 
and Roark (2010), Cohen, Safta, and Collins (2013). For a discussion of the optimality 
of the grammar constants in Valiant's algorithm, see for example Abboud, Backurs, and 
Williams (2015). 


2. Background and Notation 

This section provides background on LCFRS, and establishes notation used in the 
remainder of fhe paper. A reference fable of nofation is also provided in Appendix A. 

For an integer n, let [n] denote the set of integers {1,..., n}. Let [n]o = [n] U {0}. For 
a set X, we denote by X'^ the set of all sequences of length 1 or more of elements from 

A. 

A span is a pair of infegers denofing left and right endpoints for a substring in a 
larger string. The endpoints are placed in the "spaces" between the symbols in a string. 
For example, the span (0,3) spans the first three symbols in the string. For a string of 
length n, the set of potential endpoints is [n]o. 

We turn now to give a succinct definition for binary LCFRS. For more details about 
LCFRS and their relationship to other grammar formalisms, see Kallmeyer (2010). A 
binary LCFRS is a tuple {C, T, TZ, p, S) such that: 

• £ is the set of nonterminal symbols in the grammar. 

• T is the set of terminal symbols in the grammar. We assume CoT =0. 

• p is a function specifying a fixed fan-ouf for each nonterminal (p: £ —5- N). 

• 7^ is a set of productions. Each production p has the form A ^ g[B,C] where 
A,B,Cg C, and g is a composition function g : (T"*)*^!^) x (7"*)‘^(C') —j. (7"*)‘^(-4)^ 
which specifies how to assemble the (f{B) + (p{C) spans of fhe righfhand side 
nonferminals into the p(A) spans of the lefthand side nonterminal. We use square 
brackets as part of the syntax for writing productions, and parentheses to denote 
the application of the function g. The function g must be linear and non-erasing, 
which means that if g is applied on a pair of tuples of strings, then each input 
string appears exactly once in the output, possibly as a substring of one of the 
strings in the output tuple. Rules may also take the form A —> g[], where g returns 
a constant tuple of one string from T. 

• S' G £ is a sfarf symbol. Without loss of generality, we assume p(S) = 1. 

The language of an LCERS G = (£, T, TZ, p, S) is defined as follows: 

• We define first the set yield (A) for every A € C: 

• Eor every A ^ G TZ, g{) G yield(A). 

• Eor every A — g[B, C] €TZ and all fuples /3 G yield(i3), 7 G yield(C), 
g(/3,7) € yield(A). 

• Nothing else is in yield (A). 

• The string language of G is L{G) = {w | (w) G yield(S)}. 
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Intuitively, the process of generating a string from an LCFRS grammar consists of 
first choosing, top-down, a production to expand each nonterminal, and then, bottom- 
up, applying the composition functions associated with each production to build the 
string. As an example, the following context-free grammar: 

S^AB 
A —y a 
B 

corresponds to the following (binary) LCFRS: 


S —?■ gi[A, B] 

5i((/3i), ( 71 )) = (/3i7i) 

(1) 

^ 52 [] 

52 0 = (a) 


B 53 O 

530 = (0 



The only derivation possible under this grammar consists of the function application 

51(520, ffaO) = {ab). 

The following notation will be used to precisely represent the linear non-erasing 
composition functions g used in a specific grammar. For each production rule that 
operates on nonterminals A, B, and C, we define variables from the set S = 
{di, • • ■, 7i, • ■ •, 7<^(c)}- Iri addition, we define variables Oi for each rule where 

i G [ip(7l)], taking values from <S+. We write an LCFRS function as: 


5((0i, • ■ •, ,(7 i,*'*,7v 2(C))) (oi,..., o^ip(^A)) 


where each ai = • • • ai^m specifies the parameter strings that are combined to form 

the ith string of the function's result tuple. For example, for the rule in Eq. 1, ai,i = di 
and 01^2 = 72 - 

We adopt the following notational shorthand for LCFRS rules in the remainder of 
the paper. We write the rule: 

A^g[B,C] 

5 ((/^l, • ■ • , ), (71, * ■ * , 7 ¥,(C))) (^1, • ■ • , ) 


as: 


A[a] ^ B[I3] CU 

where a consists of a tuple of strings from the alphabet {di, ■ • ■, dv(s), 7i, • ■ •, lv{C )} • 
this notation, /? is always the tuple (di,. ■., li,p{B)), and 7 is always (71, ..., jip(c))- 
include d and 7 in the rule notation merely to remind the reader of the meaning of the 
symbols in a. 

For example, with context-free grammars, rules have the form: 

A[(di7i)]^B[(di)]C'[(7i)] 

indicating that B and C each have one span, and are concatenated in order to form A. 
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C 

B 

A 


72 


P2 


72 P2 


7i 


/3i 


/3i7i 


Figure 1 

An example of a combination of spans for Tree Adjoining Grammars (TAG) for the adjunction 
operation in terms of binary LCFRS. The rule in Eq. 2 specifies how two nonterminals B and C 
are combined together into a nonterminal A. 


A binary tree-adjoining grammar can also be represented as a binary LCFRS (Vijay- 
Shanker and Weir 1994). Figure 1 demonstrates how the adjunction operation is done 
with binary LCFRS. Each gray block denotes a span, and the adjunction operator takes 
the first span of nonterminal B and concatenates it to the first span of nonterminal C 
(to get the first span of A), and then takes the second span of C and concatenates it with 
the second span of B (to get the second span of A). For tree-adjoining grammars, rules 
have the form: 


A[(bi71,72/32)] ^ B[(/3 i,/32)]C[(7i,72)] (2) 

The fan-out of a nonterminal is the number of spans in the input sentence that it 
covers. The fan-out of CFG rules is one, and the fan-out of TAG rules is two. The fan¬ 
out of the grammar, /, is the maximum fan-out of its nonterminals: 

/ = maxp(A). (3) 

We sometimes refer to the skeleton of a grammar rule A[a\ — B[(3] Gjy], which 
is just the context-free rule A ^ B C, omitting the variables. In that context, a logical 
statement such as A ^ BC &TZ is true if there is any rule A[a] —)■ B[I3] €[7] &Tl with 
some a, (3, and 7. 

For our parsing algorithm, we assume that the grammar is in a normal form such 
that the variables /3i, ..., appear in order in a, that is, that the spans of B are 
not re-ordered by the rule, and similarly we assume that 71,..., 7(^(c) appear in order. 
If this is not the case in some rule, the grammar can be transformed by introducing a 
new nonterminal for each permutation of a nonterminal that can be produced by the 
grammar. We further assume that aij = f3i, that is, that the first span of A begins with 
material produced by B rather than by C. If this not the case for some rule, B and C can 
be exchanged to satisfy this condition. 

We refer to an LGFRS rule A ^ B C as single-initial if the leftmost endpoint of 
C is internal to a span of A, and dual-initial if the leftmost endpoint of C is the 
beginning of a span of A. Our algorithm will require the input LGFRS to be in single¬ 
initial form, meaning that all rules are single-initial. We note that grammars for common 
formalisms including TAG and synchronous context-free grammar (SGFG) are in this 
form. If a grammar is not in single-initial form, dual-initial rules can converted to 
single-initial form by adding a empty span to B which combines with the first spans 
of C immediately to its left, as shown in Figure 2. Specifically, for each dual-initial rule 
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C 


C 


B 


B' 


A 


A 


Figure 2 

Conversion of a dual-initial rule to a single-initial rule. 


A ^ B C, A the first span of C appears between spans i and i + 1 of B, create a new 
nonterminal B' with ^p{B') = <p(i?) -I- 1, and add a rule B' —t B, where B' produces B 
along with a span of length zero between spans i and i + 1 of B. We then replace the 
rule A —)■ B C with A ^ B' C, where the new span of B' combines with C immediately 
to the left of C"s first span. Because the new nonterminal B' has fan-out one greater than 
B, this grammar transformation can increase a grammar's fan-out by at most one. 

By limiting ourselves to binary LCFRS grammars, we do not necessarily restrict the 
power of our results. Any LCFRS with arbitrary rank (i.e. with an arbitrary number 
of nonterminals in the right-hand side) can be converted to a binary LCFRS (with 
potentially a larger fan-out). See discussion in §7.6. 


Example 1 

Consider the phenomenon of cross-serial dependencies that exists in certain languages. 
It has been used in the past (Shieber 1985) to argue that Swiss-German is not context- 
free. One can show that there is a homomorphism between Swiss-German and the 
alphabet {a, b, c, d} such that the image of the homomorphism intersected with the 
regular language a*b*c*d* gives the language L = {d^b^c^dA \ m,n > 1}. Since L is not 
context-free, this implies that Swiss-German is not context-free, because context-free 
languages are closed under intersection with regular languages. 

Tree-adjoining grammars, on the other hand, are mildly context-sensitive for¬ 
malisms that can handle such cross-serial dependencies in languages (where the as 
are aligned with cs and the bs are aligned with the ds). For example, a tree-adjoining 
grammar for generating L would include the following initial and auxiliary trees (nodes 
marked by * are nodes where adjunction is not allowed): 


Initial 

S 

A* 

Auxiliary 

A* 

B* 

A 

B 

a A 

/\ 
b B 

e 

A* 

A 

A* c 

A 

B* d 
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This TAG corresponds to the following LCFRS: 


S ^ gi [A] 

3i((/3i,/ 32)) = (di/32) 

^ 34 

34((/3i./ 32)) = (dl; d2) 

A —>■ 32 [A] 

32((dl,d2)) = (adi, cd2) 

B —>■ 35 [f?] 

35((dl,d2)) = {bPi, dd2) 

A 33[] 

330 = {£, s) 

B gei] 

36 0 = (£, 0 


Here we have one unary LCFRS rule for the initial tree, one unary rule for each ad- 
juncfion tree, and one null-ary rule for each nonferminal producing a tuple of empfy 
strings in order to represent TAG tree nodes at which no adjunction occurs. The LCFRS 
given above does not satisfy our normal form requiring each rule fo have eifher fwo 
nonterminals on the righthand side with no terminals in the composition function, or 
zero nonterminals with a composition function returning fixed strings of ferminals. 
However, it can be converted to such a form through a process analogous to converting 
a CFG to Chomsky Normal Form. For adjunction trees, the two strings returned by the 
composition function correspond the the material to the left and right of fhe foot node. 
The composition function merges terminals at the leaves of fhe adjuncfion tree with 
material produced by internal nodes of fhe free at which adjunction may occur. 

In general, binary LCFRS are more expressive than TAGs because they can have 
nonterminals with fan-out greater than two, and because they can interleave the argu¬ 
ments of the composition function in any order. 

3. A Sketch of the Algorithm 

Our algorithm for LCFRS sfring recognition is inspired by fhe algorifhm of Valianf 
(1975). If infroduces a few importanf novelties that make it possible to use matrix 
multiplication for the goal of LCFRS recognifion. 

The algorifhm relies on the observation that it is possible to construct a matrix T 
with a specific non-associative multiplication and addition operator such that multiply¬ 
ing T by itself k fimes on the left or on the right yields /c-step derivations for a given 
string. The row and column indices of fhe matrix together assemble a set of spans in 
fhe sfring (the fan-out of fhe grammar defermines fhe number of spans). Each cell in 
the matrix keeps track of fhe nonterminals that can dominate these spans. Therefore, 
computing the transitive closure of fhis matrix yields in each matrix cell the set of 
nonterminals that can dominate the assembled indices' spans for fhe specific string at 
hand. 

There are several key differences between Valiant's algorithm and our algorithm. 
Valiant's algorithm has a rather simple matrix indexing scheme for fhe matrix: the rows 
correspond to the left endpoints of a span and the columns correspond to its right 
endpoints. Our matrix indexing scheme can mix both left endpoints and right endpoints 
at either the rows or the columns. This is necessary because with LCFRS, spans for the 
right-hand side of an LCFRS rule can combine in various ways info a new sef of spans 
for fhe leff-hand side. 

In addition, our indexing scheme is "over-complefe." This means thaf different cells 
in the matrix T (or its matrix powers) are equivalent and should consist of fhe same 
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i j k i m n 

1 2 4 5 7 8 


[C,j,kJ,m\ 

[A,i,k,e,n\ 


Figure 3 

A demonstration of a parsing step for the combination of spans in Figure 1. During parsing, the 
endpoints of each span are instantiated with indices into the string. The variables for these 
indices shown on the left correspond to the logical induction rule on the right. The specific 
choice of indices shown at the bottom is used in our matrix multiplication example in §3. 


nonterminals. The reason we need such an over-complete scheme is again because of 
the possible ways spans of a righf-hand side can combine in an LCFRS. To address this 
over-completeness, we introduce into the multiplication operator a "copy operation" 
that copies nonterminals between cells in order to maintain the same set of nonterminals 
in equivalent cells. 

To give a preliminary example, consider the tree-adjoining grammar rule shown 
in Figure 1. We consider an application of the rule with the endpoints of each span 
instantiated as shown in Figure 3. With our algorithm, this operation will translate 
into the following sequence of matrix transformations. We will start with the following 
matrices, Ti and T 2 : 


Ti 


(2,7) 


T2 


(4,5) 



/ 


\ 


/ 


\ 

(1,8) 












(2,7) 


{...,C,...} 



1 


) 


1 


/ 


For Ti, for example, the fact that B appears for the pair of addresses (1,8) (for row) 
and (2,7) for column denotes that B spans the constituents (1,2) and (7,8) in the string 
(this is assumed to be true - in practice, it is the result of a previous step of matrix 
multiplication). Similarly, with T 2 , C spans the constituents (2,4) and (5,7). 

Note that (2,7) are the two positions in the string where B and C meet, and that 
because B and C share these two endpoints, they can combine to form A. In the matrix 
representation, (2,7) appears as the column address of B and as the row address of C, 
meaning that B and C appear in cells that are combined during matrix multiplication. 
The result of multiplying Ti by T 2 is the following: 
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TiT2 


(4,5) 



/ 


\ 

(1,8) 





1 


/ 


Now A appears in the cell that corresponds to the spans (1,4) and (5,8). This is the 
result of merging the spans (1,2) with (2,4) (left span of B and left span of C) into (1,4) 
and the merging of the spans (5, 7) and (7,8) (right span of C and right span of B) into 
(5,8). Finally, an additional copying operation will lead to the following matrix: 


T3 


(5,8) 


/ 

\ 

(1,4) 1 




V ; 

Here, we copy the nonterminal A from the address with the row (1,8) and column 
(4,5) into the address with the row (1,4) and column (5,8). Both of these addresses cor¬ 
respond to the same spans (1,4) and (5,8). Note that matrix row and column addresses 
can mix both starting points of spans and ending points of spans. 

4. A Matrix Multiplication Algorithm for LCFRS 

We turn next to give a description of the algorithm. Our description is constructed as 
follows: 

• In §4.1 we describe the basic matrix structure which is used for LCFRS recognition. 
This construction depends on a parameter d, the contact rank, which is a function 
of the underlying LCFRS grammar we parse with. We also describe how to create 
a seed matrix, for which we need to compute the transitive closure. 

• In §4.2 we define the multiplication operator between cells of the matrices we use. 
This multiplication operator is distributive, but not associative, and as such, we 
use Valiant's specialized transitive closure algorithm to compute transitive closure 
of the seed matrix given a string. 

• In §4.3 we define the contact rank parameter d. The smaller d is, the more efficient 
it is to parse with the specific grammar. 

• In §4.4 we define when a binary LCFRS is “balanced." This is an end case that in¬ 
creases the final complexity of our algorithm by a factor of 0{n). Nevertheless, it is 
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an important end case that appears in applications, such as inversion transduction 
grammars. 

• In §4.5 we tie things up, and show that computing the transitive closure of the 
seed matrix we define in §4.1 yields a recognition algorithm for LCFRS. 

4.1 Matrix Structure 

The algorithm will seek to compute the transitive closure of a seed mafrix T{d), where 
d is a consfant defermined by fhe grammar (see §4.3). The matrix rows and columns are 
indexed by the set N{d) defined as: 


d 

(4) 

d’=l 

where n denofes fhe lengfh of fhe senfence, and the exponent d' denotes a repeated 
Cartesian product. Thus each element of N{d) is a sequence of indices info the string, 
where each index is annotated with a bit (an element of the set {0,1}) indicating whether 
it is marked or unmarked. Marked indices will be used in the copy operator defined 
lafer. Indices are unmarked unless specified as marked: we use x fo denote a marked 
index (x, 1) with x £ [n]o. 

In the following, if will be safe to assume sequences from N (d) are monotonically 
increasing in fheir indices. For an f s (d), we overload nofation, and often refer to the 
set of all elements in the first coordinate of each elemenf in fhe sequence (ignoring fhe 
additional bits). As such, 

• The set f U j is defined for j £ N (d). 

• If we stafe that f is in N{d) and includes a set of endpoinfs, if means fhat f is a se¬ 
quence of these integers (ordered lexicographically) with the bit part determined 
as explained in the context (for example, all unmarked). 

• The quanfity |f | denotes the length of the sequence. 

• The quantity min i denotes the smallest index among the first coordinates of all 
elements in the sequence i (ignoring the additional bits). 

We emphasize that the variables i, j, and k are mostly elements in N (d) as over¬ 
loaded above, not integers, throughout this paper; we choose the symbols i, j, and k 
by analogy to the variables in the CKY parsing algorithm, and also because we use the 
sequences as addresses for matrix rows and columns. For i,j £ N(d), we define m{i, j) 

tobefhesef of/' = i|fU j| pairs {(£ 1 ,^ 2 ), (4,^4), 4/')} such fhat 4 < 4+i 

for k £ [2/' — 1] and (ffc, 0) £ iU j for k £ [2/']. This means thaf m{i,j) takes as input 
the two sequences in matrix indices, merges them, sorts them, then divides this sorted 
list into a set of /' consecutive pairs. Whenever min j < mini, m{i,j) is undefined. The 
inferprefation of this is that f'l should always belong to i and not j. See more details in 
§4.2. In addition, if any element of i or j is marked, is undefined. 

We define an order < on elements i and j of N (d) by first sorting the sequences 
i and j and then comparing i and j lexicographically (ignoring the bits). This ensures 
that z < j if mini < minj. We assume fhat the rows and columns of our matrices are 
arranged in this order. For the rest of the discussion, we assume that d is a constant, and 
refer fo T{d) as T and N (d) as N. 
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We also define the set of triples M as the following Cartesian product: 

M = (£ U {[i, [h, [h, !;□, D}) X N X N, 

where ci, Lh, Lh, a, and □ are six special pre-defined symbols.^ Each cell in T is a 
set such that c M. 

The intuition behind matrices of the type of T (meaning T and, as we see later, 
products of T with itself, or its transitive closure) is that each cell indexed by (z, j) in 
such a matrix consists of all nonterminals that can be generated by the grammar when 
parsing a sentence such that these nonterminals span the constituents (whenever 

is defined). Our normal form for LCFRS ensures that spans of a nonterminal are 
never re-ordered, meaning that it is not necessary to retain information about which 
indices demarcate which components of the nonterminal, because one can sort the 
indices and take the first two indices as delimiting the first span, the second two indices 
as delimiting the second span, and so on. The two additional N elements in each triplet 
in a cell are actually just copies of the row and column indices of that cell. As such, 
they are identical for all triplets in that cell. The additional i, Lh, Lh, i, a, □ symbols are 
symbols that indicate to the matrix multiplication operator that a "copying operation" 
should happen between equivalent cells (§4.2). 

Figure 4 gives an algorithm to seed the initial matrix T. Entries added in step 2 
of the algorithm correspond to entries in the LCFRS parsing chart that can be derived 
immediately from terminals in the string. Entries added in step 3 of the algorithm do 
not depend on the input string or input grammar, but rather initialize elements used 
in the copy operation described in detail in §4.2. Because the algorithm only initializes 
entries with z < j, the matrix T is guaranteed to be upper triangular, a fact which we 
will take advantage of in §4.2. 

4.1.1 Configurations. Our matrix representation requires that a nonterminal appears 
in more than one equivalent cell in the matrix, and the specific set of cells required 
depends on the specific patterns in which spans are combined in the LCFRS grammar. 
We now present a precise description of these cells by defining the configuration of a 
nonterminal in a rule. The concept of a configuration is designed to represent which 
endpoints of spans of the rule's righthand side (r.h.s.) nonterminals B and C meet 
one another to produce larger spans, and which endpoints, on the other hand, become 
endpoints of spans of the lefthand side (l.h.s.) nonterminal A. 

For each of the three nonterminals involved in a rule, the configuration is the set of 
endpoints in the row address of the nonterminal's matrix cell. To make this precise, for 
a nonterminal B with fan-out (p{B), we number the endpoints of spans with integers in 
the range 1 to 2(p{B). In a rule A[a\ —>■ i?[/3] C'[ 7 ], the configuration of B is the subset 
of [2(p{B)] of endpoints of B that do not combine with endpoints of C in order to 
form a single span of A. The endpoints will form the row address for B. Formally, let 
(3 = (/3i,... and let a = (ai,i • • Then the set 

of non-combination endpoints of B is defined as: 

config 2 (r) = {2z | Pi = aj^n, for some j} U {2z — 1 | /3j = for some j} 


2 The symbols will be used for "copying commands:" (1) "from row" (A); (2) "from column" (O); (3) "to 
row" (i); (4) "to column" ([b); (5) "unmark row" (□); (6) "unmark column" (D). 
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Inputs: An LCFRS grammar as defined in §2 and a sentence wi ■ ■ -Wn- 

Outputs: A seed matrix T with rows and columns indexed by N, such that each cell in 

T is a subset of M 

Algorithm: 

1. Set Tij = 0 for all f, j G N. 

2 . For each i,j€N, for each nonterminal A G C, set U {(A, i, j)} if 

f^ 2 )) (^3 5 ) (^ 2 /-i)^ 2 /)} arid there is a rule in the grammar A —)■ 

90, 5() = (ai, ■ • ■, where a* = and for each i and j, aij = 

rr^ 2 i-i+i (i-S- “ij is iias (^ 2 i-i + i)th word in the sentence). 

3. For each i,j G N such that i < j 

a. Tij ■(— Tij U {([F, i, j)} if all indices in i are unmarked and j = insert(f, x) for 
some X 

b. Tij Tij U {(i, f, j)} if all indices in j are unmarked and i = remove(j, x) for 
some X 

c. Ty ■(—Ty U {( 1 = 1,15 j)} if all indices in i are unmarked and i = 
insert (remove(j, x), x) for some x 

d. Tij Tij U {(i, i, j)} if all indices in j are unmarked and i = insert(j, x) for 
some X 

e. Tij ■(— Tij U {([h, i, j)} if all indices in i are unmarked and j = remove(i, x) for 
some X 

f. Ty ■(—Ty U {(□, h j)} if all indices in j are unmarked and j = 
insert (remove(i, x), x) for some x 


Figure 4 

An algorithm for computing the seed matrix T. The function remove(w, x) takes a sequence of 
integers v and removes x from it, if it is in there. The function insert(i;, x) takes a sequence of 
integers and adds x to it. 


where the first set defines right ends of spans of B that are right ends of some span of 
A, and the second set defines left ends of spans of B that are left ends of some span of 
A. For example, given that CFG rules have the from 

r = A[(/3 i7i)] ^ B[{f3i)] G[( 7 i)] 

the configuration config 2 (r) is {1} because, of B's two endpoints, only the first is also an 
endpoint of A. For the TAG rule t shown in Figure 1, config 2 (f) = {1,4} because, of B's 
four endpoints, the first and fourth are also endpoints of A. 

For the second r.h.s. nonterminal of a rule r, the configuration consists of the set of 
endpoints in the row address for C, which are the endpoints that do combine with B: 

config 3 (r) = {2i \ 7^ = aj^k for some 1 < fc < Uj} U {2i — 1 | 7^ = 0^7, for some 1 < k < rij} 

where the first set defines right ends of spans of C that are internal to some span of 
A, and the second set defines lefts ends of spans of C that are internal to some span 
of A. For example, any CFG rule r has configuration, config 3 (r) = {1}, because the first 
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endpoint of C is internal to A. For the TAG rule t shown in Figure 1, config 3 (t) = {1,4} 
because, of G's four endpoints, the first and fourth are internal A. 

For the l.h.s. nonterminal A of the rule, matrix multiplication will produce an entry 
in the matrix cell where the row address corresponds to the endpoints from B, and the 
column address corresponds to the endpoints from C. To capture this partition of the 
endpoints of A, we define 

config 3 (r) = {2i | = Pj for some jj U {2i — 1 | = jSj for some j}, 

where the first set defines right ends of spans of A that are formed from B, and the 
second set defines left ends of spans of A that are formed from B. For example, any 
CFG rule r has configuration, config^(r) = {1}, because only the first endpoint of A is 
derived from B. For the TAG rule t shown in Figure 1, config^ (t) = {1,4} because, of A's 
four endpoints, the first and fourth are derived from B. 

4.2 Definition of Multiplication Operator 

We need to define a multiplication operator 0 between a pair of elements R, S C M. 
Such a multiplication operator induces multiplication between matrices of the t 5 rpe of 
T, just by defining for two such matrices, Ti and T 2 , a new matrix of the same size 
Ti 0 T 2 such that: 


[Ti 0 T 2 ]ij= [J {[Ti]ik 0 [T 2 ]kj), (5) 

fesAf 

We also use the U symbol to denote coordinate-wise union of cells in the matrices it 
operates on. 

The operator 0 we define is not associative, but it is distributive over U. This means 
that for R, 81,82 0 M it holds that: 

i?(g) (S'! U 82 ) = {R0 81 ) U (i?(g 82 ). 

In addition, whenever i? = 0, then for any 8 , R 0 8 = 8 0 R = ^. This property 
maintains the upper-triangularity of the transitive closure of T. 

Figure 5 gives the algorithm for multiplying two elements of the matrix. The algo¬ 
rithm is composed of two components. The first component (step 2 in Figure 5) adds 
nonterminals, for example. A, to cell (i, j), if there is some B and C in {i,k) and (fc, j), 
respectively, such that there exists a rule A ^ B C and the span endpoints denoted by k 
are the points where the rule specifies that spans of B and C should meet. 

In order to make this first component valid, we have to make sure that k can 
indeed serve as a concatenation point for {i,j). Step 2 verifies this using the concept 
of configurations defined above. To apply a rule r : A[a] —?■ B[P] Gjy], we must have an 
entry for {B, i, k) in cell (1, k), where i is a set of indices corresponding to the endpoints 
of B selected by config 2 (r) and fc is a set of indices corresponding to the endpoints of 
B selected by [2if{B)] \ config 2 (r). This condition is enforced by step 2c of Figure 5. 
Similarly, we must have an entry for {C, k,j) in cell {k,j), where fc is a set of indices 
corresponding to the endpoints of C selected by config 3 (r) and j is a set of indices 
corresponding to the endpoints of C selected by [2(p(G)] \ config 3 (r). This is enforced 
by step 2d. Finally, the spans defined by B and C must not overlap in the string. To 
guarantee that the spans do not overlap, we sort the endpoints of A and check that each 
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Inputs: A pair of elements R,Sc M. 

Outputs: A new subset of M, denoted by {R® S). 

Algorithm: 

1. {R®S) = %. 

2. For each pair of elements r = {B, i, k) and s = (C, k, j) where r G R and s G S 
add an element {A, i,j) to (i? 0 S') if: 

a. There is a binary rule in the LCFRS grammar r = A[a] —>■ B[/3] C['y]. 

b. All indices in i, j, and k are unmarked. 

c. Configuration of B: assume m{i, k) = , (4v(b)-u4v(b))}- 

We require that i = {£q \ q G config 2 (r)} 

d. Configuration of C: assume m{k,j) = {(^ 1 ,^ 2 )) ■ • ■) (^ 2 <, 3 (C)-i> ^2<p(C))}' 

We require that k = {i'q \ q G config 3 (r)} 

e. Configuration of A: assume > (.^ 2 v(a)-v^ 2 v{A))}- 

We require that i = {£" | q G config^(r)} 

3. For each pair of elements r = (A, i, k) and s = (IF, k, j), add the element (A, i, j) 
to (i? (8> S) if: i contains only unmarked indices, and x G j and x G i for some x. 

4. For each pair of elements r = (i, i, k) and s = (A, fc, j), add the element (A, i, j) 
to (i? (8> S) if: X G j and x ^ i and for some x. 

5. For each pair of elements r = (A, i, k) and s = (□, k,j), add the element (A, i,j) to 
{R 0 S) if: |t U j| = 2(p(A). 

6. For each pair of elements r = (i, i, k) and s = (A, fc, j), add the element (A, i, j) 
to (i? (8> S) if: j contains only unmarked indices, and x G i and x G j for some x. 

7. For each pair of elements r = (A, i, k) and s = (LF, k, j), add the element (A, i, j) 
to (i? (g) S) if: a; € t and x ^ j for some x. 

8. For each pair of elements r = ( 1 = 1 , i, k) and s = (A, k,j), add the element (A, i, j) 
to (i?® S) if: \i U j| = 2ip(A). 


Figure 5 

An algorithm for the product of two matrix elements. 


position in the sorted list is derived from either B or C as required by the configuration 
of A in r. This check is performed in step 2e of Figure 5. 

Given that T is initialized to be upper-triangular, the properties of matrix multipli¬ 
cation guarantee that all matrix powers of T are upper-triangular. We now proceed to 
show that upper-triangular matrices are sufficient in terms of the grammar. In particular, 
we need to show the following lemma: 

Lemma 1 

For each application of a single-initial rule A —> i? G, it is possible to create an entry for 
A by multiplying two upper-triangular matrices Ti and T 2 , where Ti contains an entry 
for B, and T 2 contains an entry for C. 

Proof 

A nonterminal B appears in a cell above the diagonal if its row address is smaller 
than its column address, which in turn occurs if the leftmost endpoint of B appears 
in the row address rather than the column address. The row address for B contains the 
endpoints of B that are also endpoints of A. Our normal form for LCFRS rules ensures 
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that the leftmost endpoint of B forms the leftmost endpoint of A. Therefore fhe leffmosf 
endpoinf of B is in B's row address, and B is above fhe diagonal. 

The row address of nonterminal C in T 2 must contain the endpoints of C thaf com¬ 
bine wifh endpoinfs of B. For single-initial rules, fhese endpoints include the leftmost 
endpoint of C, guaranfeeing that C appears above the diagonal. 

Because each instance of A can be produced by combining elements of Ti and T 2 
fhat are above fhe diagonal, each insfance of A can be produced by mulfiplying fwo 
upper-friangular matrices. ■ 

4.2.1 Copy Operations. The first component of the algorithm is sound, but not complete. 
If we were fo use just this component in the algorithm, then we would get in each cell 
{i,j) of the transitive closure of T a subset of the possible nonterminals that can span 
The reason this happens is that our addressing scheme is "over-complete." This 
means that any pair of addresses (f, j) and {k, i) are equivalent if m{i,j) = m{k, £). 

We need fo ensure thaf the transitive closure, using (g), propagates, or copies, non¬ 
terminals from one cell to its equivalents. This is done by the second component of the 
algorithm, in steps 3-6. The algorithm does this kind of copying by using a set of six 
special "copy" symbols, 0 , 1 = 1 }. These symbols copy nonferminals from one 

cell fo the other in multiple stages. 

Suppose that we need to copy a nonterminal from cell {i,j) to cell (fc,f), where 
’m-ihj) = m{k,i), indicating that the two cells describe the same set of indices in the 
input string. We must move the indices in f n £ from fhe row address fo the column 
address, and we must move the indices in j Dk from the column address to the row 
address. We will move one index at a time, adding nonterminals to intermediate cells 
along the way. 

We now illustrate how our operations move a single index from a row address to 
a column address (moving from column fo row is similar). Lef x indicate the index 
we wish to move, meaning that we wish to copy a nonterminal in cell {i,j) to cell 
(remove(z,x), insert(j, x)). Because we want our overall parsing algorithm to take ad¬ 
vantage of fasf matrix multiplication, we accomplish the copy operations through a 
sequence of fhree mafrix multiplications, as shown in Figure 6. The first multiplication 
involves the nonterminal A in cell (f,j) in the left matrix, and a LF symbol in cell 
(j, insert)}, x)) in the right matrix, resulting in a matrix with nonterminal A in cell 
{i, insert)}, x)). This intermediate result is redundant in the sense that index x appears in 
the row and index x appears in the column address. To remove x from the row address, 
we multiply on the left with a matrix containing the symbol i in cell (remove)^, x), i), 
resulting in a matrix with nonterminal A in cell (remove)^, x), insert)}, x)). Finally, we 
multiply by a third matrix to replace the marked index x with the unmarked index x. 
This is done by multiplying on the right with a matrix containing the symbol D in cell 
(insert)}, x), insert)}, x)). 

The key idea behind the above three-step process is to copy elements from one cell 
fo anofher through intermediate cells. In matrix multiplication, only cells that share a 
row or a column index actually interact when doing multiplication. Therefore, in order 
fo copy a nonterminal from (f,}) fo another cell which represents the same set of spans, 
we have to copy it through cells such as {i, insert)}, x)) that share the index i with (f, }). 

In order to guarantee that our operations copy nonterminals only into cells with 
equivalent addresses, the seed matrix contains the special symbol LF only in cells (}, k) 
such that k = insert)}, x) for some x. When [F in cell (}, k) combines wifh a nonferminal 
A in cell (f, }), the result contains A only if x G f, guaranfeeing fhat the index added to 
the column address was originally present in the row address. In addition, the condition 
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Figure 6 

An example of moving an index from the row address to the column address. Nonterminal B in 
Ti is copied from cell (1, 8 ), (2, 7) to cell (1), (2, 7, 8 ) through three matrix multiplications. First, 
multiplying by T 2 on the right yields T 1 T 2 , shown in the right of the second row. Multiplying 
this matrix by T 3 on the left yields T 1 T 2 T 3 . Finally, multiplying this matrix by T 4 on the right 
yields TiT 2 T 3 Tn, shown in the bottom row. 
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that i contains only unmarked indices (in the multiplication operator) and that the 
condition j contains only unmarked indices (in the initialization of the seed matrix) 
guarantee that only one index is marked in the address of any non-empty matrix cell. 

Similar conditions apply to the operation. The seed matrix contains only in 
cells (i, k) such that i = remove(fc, x) for some x, guaranteeing that the operation only 
removes one index at a time. Furthermore, when in cell (z, k) combines with a 
nonterminal A in cell {k,j), the result contains A only if x € j. This guarantees that 
the new entry includes all the original indices, meaning that any index we remove from 
the row address is still present as a marked index in the column address. 

The □ operator removes the mark on index x in the column address, completing 
the entire copying process. The condition |z U j | = ip{A) ensures that the removal of the 
mark from x does not take place until after x has been removed from the row address. 

Taken together, these conditions ensure that after a sequence of one Lh, one ci, and 
one D, A is copied into all cells having the form (remove(z, x), insert(j, x)) for some x. 

To move an index from the column address to the row address, we use one [i 
operation followed by one Ih operation and one i=i operation. The conditions on these 
three special symbols are analogous to the conditions on [F, i, and D outlined above, 
and ensure that we copy from cell (z, j) to cells of the form (insert(z, x), remove(j, x)) for 
some X. 

We now show that matrix powers of the upper-triangular seed matrix T copy 
nonterminals between all equivalent cells above the diagonal. 

Lemma 2 

Let (z,j) and {k,t) be unmarked matrix addresses, in a seed matrix T indexed by 
row and column addresses from N{d) where d > min{|z|, |_)|} and d > min{|fc|, \^\}. 
Assume that mini = minfc and either k = remove(z,a;) and i = insert(j,a;) for some x, 
or k = insert(z, x) and £ = remove(j, x) for some x. If A appears in cell (z, j) of T^'^\ then 
A appears in cell {k,tj of t("+ 3). Furthermore, the copy operations do not introduce 
nonterminals into any other cells with unmarked addresses. 

Proof 

The condition on d guarantees that we can form row and column addresses long 
enough to hold the redundant representations with one address shared between row 
and column. This condition is only relevant in the case where z, j, k, and I are all of 
the same length; in this case we need to construct temporary indices with length one 
greater, as in the example in Figure 6. 

A can be added to cell (/c, through a sequence of three matrix multiplications 
by combining with symbols [F, and D or with i, [F, and i=i. Because is upper 
triangular, min z = min z U j, meaning that A’s leftmost index is in its row address. The 
condition min z = min k implies that we are not moving this leftmost index from row to 
column. The addresses of the three copy symbols required are all formed by adding or 
removing a; or i: to the row and column addresses (z, j); because the leftmost index of z 
is not modified, the copy symbols that are required are all above the diagonal, and are 
present in the seed matrix T. Therefore, A appears in cell (fc, £) of 

To see that nonterminals are not introduced into any other cells, observe that i=i and 
D are the only symbols that introduce nonterminals into unmarked addresses. They can 
only apply when a marked index is present, and when the total number indices is 2ip{A). 
This can only occur after either IF has introduced a marked index and ci removed the 
corresponding unmarked index, or ci has introduced a marked index and O removed 
the corresponding unmarked index ■ 
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Putting together sequences of these operations to move indices, we get the follow¬ 
ing lemma: 

Lemma 3 

Let {i,j) and {k,£) be matrix addresses such that m{i,j) = m{k,£), in a seed matrix 
T indexed by row and column addresses from N{d) where d > min{|i|, |j|} and d > 
min{|A:|, |f|}. Then, for any nonterminal A in cell {i,j) in T^”^, A will also appear in cell 
(fc, £) of the power matrix 

Proof 

Nonterminal A can be copied through a series of intermediate cells by moving one index 
at a time from i to £, and from j to k. We begin by moving indices from either the row 
address i to the column address if |i| > |j|, or from the column address j to the row 
address otherwise. We must move up to d indices from row to column, and d indices 
from column to row. Each move takes three matrix multiplications, for a total of 6 d 
matrix multiplications. ■ 

4.3 Determining the Contact Rank 

The dimensions of the matrix T (and its transitive closure) are |iV| x \N\. The set N is 
of size 0{riA), where d is a function of the grammar. When a given pair of cells in two 
matrices of the type of T are multiplied, we are essentially combining endpoints from 
the first multiplicand column address with endpoints from the second multiplicand row 
address. As such, we have to ensure that d allows us to generate all possible sequences 
of endpoints that could potentially combine with a given fixed LCFRS. 

We refer to the endpoints at which a rule's r.h.s. nonterminals meet as combining 
points. For example, in the simple case of a CFG with a rule S — NP VP, there is one 
combining point where NP and VP meet. For the TAG rule shown in Figure 1, there 
are two combining points where nonterminals B and C meet. For each rule r in the 
LGFRS grammar, we must be able to access the combining points as row and column 
addresses in order to apply the rule with matrix multiplication. Thus, d must be at least 
the maximum number of combining points of any rule in the grammar. The number of 
combining points 5{r) for a rule r can be computed by comparing the number of spans 
on the l.h.s. and r.h.s. of the rule: 

<i(A[a] ^ B[P\ CU) = p(G) + p(i?) - p(A). (6) 

Note that 6 {r) depends only on the skeleton of r (see §2), and therefore it can be 
denoted hy 5{A ^ B C)? 

For each nonterminal on the r.h.s. of the rule, the address of its matrix cell consists of 
the combination points in one dimension (either row or column), and the other points in 
the other dimension of the matrix. For r.h.s. nonterminal B in rule A^ B C, the number 
of non-combination endpoints is: 

2ip{B) - 5{A ^ B C). 


3 To see that Eq. 6 is true, consider that if we take ip{B) + </j(C) variables from the spans of the r.h.s. and 
try to combine them together to ^^(.4) sequences per span of the l.h.s., we will get ip(i?) -I- ^p{C) — ip(A) 
points where variables "touch." If <p(A) = 1, then this is clearly true. For ip(A) > 1, consider that for each 
span, we "lose" one contact point. 


18 



Cohen and Gildea 


Parsing LCFRS with Fast Matrix Multiplication 


Thus, taking the maximum size over all addresses in the grammar, the largest 
addresses needed are of length: 


r 5{A^BC), ] 

d = max max < 2(p{B) — S(A —t B C), > . 

\2^{C)-5{A-, BC) \ 

We call this number the contact rank of the grammar. As examples, the contact rank 
of a CFG is one, while the contact rank of a TAG is two. A simple algebraic manipulation 
shows that the contact rank can be expressed as follows: 

r p(A)+p(B)-p(G), ] 
d= max max < Lp{A) — ip{B) + tp{C), >. 

\-p(A) + p(i?) + p(C)J 

We require our grammars to be in single-initial form, as described in §2. Because 
the process of converting an LCFRS grammar to single-initial form increases its fan-out 
by at most one, the contact rank is also increased by at most one. 


4.4 Balanced Grammars 


We define the configuration set of a nonterminal A to the the set of all configurations 
(§4.1.1) in which A appears in a grammar rule, including both appearances in the r.h.s. 
and as the l.h.s. 


config(a) = [J {config;i(r)} U u {config 2 (r)} U |J {config 3 (r)} 


k r:lhs(r)=A 


t r:rhsl(r)=A 


ir:rhs2{r)—A 


For example, in a CFG, the configuration set of any nonterminal is {{1}}, because, as 
shown in §4.1.1, nonterminals are always used in the unique configuration {1}. For 
TAG, the configuration set of any nonterminal is {{1,4}} because, as in CFG, nontermi¬ 
nals are always used in the same configuration. 

A configuration c of nonterminal B is balanced if |c| = (p{B). This means that the 
number of contact points and non-contact points are the same. 

The contact rank d defined in the previous section is the maximum size of any 
configuration of any nonterminal in any rule. For a given nonterminal B, if (f{B) < d, 
then we can copy entries between equivalent cells. To see this, suppose that we are 
moving from cell {i,j) to (fc, t) where the length of i is greater than the length of j. As 
long as we move the first index from row to column, rather than from column to row, 
the intermediate results will require addresses no longer than the length of i. 

However, if (p{B) = d, then every configuration in which B appears is balanced: 

Vc € config(i?) |c| = (p{B) 

If <f{B) = d and B appears in more than one configuration, that is, |config(i?)| > 1, it is 
impossible to copy entries for B between the cells using a matrix of size {2n)‘^. This is 
because we cannot move indices from row to column or from column to row without 
creating an intermediate row or column address of length greater than d as a result of 
the first [K or ci operation. 


19 



Volume XX, Number xx 


We define a balanced grammar to be a grammar containing a nonterminal B such 
that ‘^{B) = d, and |config(_B)| > 1. As examples, a CFG is not balanced because, while, 
for each nonterminal B, ip{B) = d = 1, the number of configurations |config(i?)| is one. 
Similarly, TAG is not balanced, because each nonterminal has only one configuration. 
Inversion Transduction Grammars (ITGs) are balanced, because, for each nonterminal 
B, ip{B) = d = 2, and nonterminals can be used in two configurations, corresponding 
to straight and inverted rules. 

The following condition will determine which of two alternative methods we use 
for the top level of our parsing algorithm. 

Condition 4.1 

Unbalanced Grammar Condition There is no nonterminal B such that '^{B) = d and 
|config(i3)| > 1. 

This condition guarantees that we can move nonterminals as necessary with matrix 
multiplication: 

Lemma 4 

Let {i,j) and {k, £) be matrix addresses such that m{i,j) = m{k, £). Under Condition 4.1, 
for any nonterminal A in cell {i, j) in T^'^\ A will also appear in cell {k, £) of the power 
matrix r("+6<i). 


Proof 

The number of A's endpoints is 2ip{A) = |i| + |j| = |A:| + \£\. If the grammar is not bal¬ 
anced, then d > ‘^{A), and therefore d > min{|i|, |j|} and d > min{|fc|, \£\}. By Lemma 3, 
A will appear in cell (fc, () of the power matrix ■ 

4.5 Computing the Transitive Closure of T 

The transitive closure of a matrix T is the result of repeated applications of the matrix 

multiplication operator described in Eq. 5. With T being the seed matrix, we define 

T+ = y(l) y 7^(2) y ^ 

where is defined recursively as: 

2 ^( 1 ) _ 2 ^ 


2-1 

rW = (J , 

i=i 

Under Condition 4.1, one can show that given an LCFRS derivation tree t 
over the input string, each node in t must appear in the transitive closure matrix 
T+. Specifically, for each node in t representing nonterminal A sparming endpoints 
{(^ 1 . 4 ), (4,4), • • •, 4v(A))}/ at each cell TA in the matrix such that m{i,j) = 

{{£i,£ 2), (£3, £4), ■ ■ ■, (£2lp(A)-l,£2^p(A))}^ contains A. This leads to the following result: 

Lemma 5 

Under Condition 4.1, the transitive closure of T is such that represents the set of 

nonterminals that are derivable for the given spans in m{i, j). 
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Proof 

The proof is by induction over the length of the LCFRS derivations. By Lemma 1, 
derivations consisting of a single rule A\a] ^[ 7 ] produce A € for i and 

j corresponding the non-combination points of B and C. For all other i and j such 

that m{i,j) = {(£i,4), (4,4), • ■ •, (4 v(A)-i,4</,(A))}, an entry is produced in by 
Lemma 4. By induction, ^*( 6 ^+ 2 ) contains entries for all LCFRS derivations of depth s, 
and T+ contains entries for all LCFRS derivations of any length. 

In the other direction, we need to show that all entries A in T+ correspond to a valid 
LCFRS derivation of nonterminal A spanning endpoints This can be shown 

by induction over the number of matrix multiplications. During each multiplication, 
entries created in the product matrix correspond either to the application of an LCFRS 
rule with l.h.s. A, or to the movement of an index between row and column address for 
a previously recognized instance of A. ■ 

The transitive closure still yields a useful result, even when Condition 4.1 does not 
hold. To show how it is useful, we need to define the "copying" operator, FI, which takes 
a matrix T' of the same type of T, and sets n(T') using the following procedure: 

1. Define e{i,j) = \ m{i',j') = m{i,j)}, i.e. the set of equivalent configura¬ 

tions to (i, j). 

2. Set [U{T%j = U A,,,. 

This means that FI takes a completion step, and copies all nonterminals between all 
equivalent addresses in T'. Note that the FI operator can be implemented such that it 
operates in time All it requires is just taking 0{n‘^) unions of sets (corresponding 

to the sets of nonterminals in the matrix cells), where each set is of size 0 ( 1 ) with respect 
to the sentence length (i.e. the size is only a function of the grammar). 

This procedure leads to a recognition algorithm for binary LCFRS that do not satisfy 
Condition 4.1 (we also assume that these binary LCFRS would not have unary cycles or 
e rules). This algorithm is given in Figure 9. It operates by iterating through transitive 
closure steps and copying steps until convergence. When we take the transitive closure 
of T, we are essentially computing a subset of the derivable nonterminals. Then, the 
copying step (with 11) propagates nonterminals through equivalent cells. Now, if we 
take the transitive closure again, and there is any way to derive new nonterminals 
because of the copying step, the resulting matrix will have at least one new nonter¬ 
minal. Otherwise, it will not change, and as such, we recognized all possible derivable 
nonterminals in each cell. 

Lemma 6 

For any single-initial LCFRS, when step 2 of the algorithm in Figure 9 converges, T is 
such that [T]ij represents the set of nonterminals that are derivable for the given spans 
in 

Proof 

Any LCFRS derivation of a nonterminal can be decomposed into a sequence of rule 
applications and copy operations, and by induction over the length of the derivation, 
all derivations will be found. Each matrix operation only produces derivable LCFRS 
nonterminals, and by induction over the number of steps of the algorithm, only deriv¬ 
able nonterminals will be found. ■ 
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Figure 7 

Reduction of transitive closure to Boolean matrix multiplication. Boolean matrix operations 
implementing the matrix multiplication example of Section 3. 


4.5.1 Reduction of Transitive Closure to Boolean Matrix Multiplication. Valiant 
showed that his algorithm for computing the multiplication of two matrices, in terms 
of a multiplication operator similar to ours, can be reduced to the problem of Boolean 
matrix multiplication. His transitive closure algorithm requires as a black box this two- 
matrix multiplication algorithm. 

We follow here a similar argument. We can use Valiant's algorithm for the compu¬ 
tation of the transitive closure, since our multiplication operator is distributive (with 
respect to U). To complete our argument, we need to show, similarly to Valiant, that the 
product of two matrices using our multiplication operator can be reduced to Boolean 
matrix multiplication. 

Consider the problem of multiplication a matrix Ti and T 2 , and say Ti®T 2 = T^. To 
reduce it to Boolean matrix multiplication, we create 2|77.| pairs of matrices, Gr and Hr, 
where r ranges over TZ. The size of Gr and Hr is N x N.lir = A[cx\ —>• i?[/3]C'[7], we set 
[Gr\ik to be 1 if the nonterminal B appears in \Ti]ik and B, i, and k meet the conditions 
of step 2c of Figure 5. Similarly, we set [Hr]kj to be 1 if the nonterminal C appears in 
\T2]kj and C, k, and j meet the conditions of step 2d. All other cells, in both Gr and Hr, 
are set to 0. Note that Gr and Hr for all r G 77 are upper triangular Boolean matrices. 

In addition, we create 2|£| pairs of matrices, Ga and Ha, where A ranges over the 
set of nonterminals C. We set [GA]ik to be 1 if the nonterminal A appears in [Ti]ik, 
regardless the conditions of step 2c of Figure 5. Similarly, we set [HA]kj to be 1 if the 
nonterminal A appears in [T 2 ]kj, regardless of the conditions of step 2d. All other cells, 
in both Ga and Ha, are set to 0. Again, Ga and Ha for all A G £ are upper triangular 
Boolean matrices. 
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Finally, we create six additional matrices, for each element in the set 
{[h, D, □}. These matrices indicate the positions in which each symbol appears 

in the seed matrix T defined in Figure 4: 

1. G^, for which [G^ij = 1 only if i,j) G T. 

2. iF[p, for which = 1 only if (IF, i,j) &T. 

3. Hq, for which [H^ij = 1 only if (□, i,j) G T. 

4. G^, for which [G^ij = 1 only if (i, i,j) G T. 

5. -ff[p, for which [H^ij = 1 only if (LF, i, j) e T. 

6 . Gen, for which [G^^ij = 1 only if (i=i, i,j) G T. 

Now, for each rule r € TZ, we compute the matrix G = GrHr- The total number of 
mafrix multiplications required is \R\, which is constant in n. Now, can be obtained 
by multiplying these matrices, and applying the conditions of Figure 5: 

1. For each A & L, for each rule r = A ^ B G, check whefher [Ir]ij = 1. If sfep 2e is 
satisfied for A, i, and j, fhen add {A, i,j) to [Ta]^. 

2. For each A G C, compute Ja = GaH^. For each (i, j), add A to [Tsjy if i G j and 
X €i for some x, and [JaIij = 1- 

3. For each A G C, compufe Ja = G^Ha- For each {i,j)> add A to [Tsjy x G j and 
X ^ i and for some x, and = 1- 

4. For each A G C, compufe Ja = GaH^]- For each add A to [Ta]^ if \i U j| = 
2lp{A), and [JA]tj = 1- 

5. For each Ag C, compute Ja = G^Ha- For each {i, j), add A to if x G i and 
X G j for some x, and [JA]ij = 1- 

6. For each A G C, compufe Ja = GaH\^- For each (z, j), add A to [Ta]^ x G i and 
X ^ j for some x, and [JA\ij = 1- 

7. For each Ag G, compufe Ja = G^Ha- For each (z, j), add A to [Ta]^ if |z U j| = 
2(p(A), and [JA]rj = 1- 


Lemma 7 

The mafrix product operation for two matrices of size {2n)‘^ x (2n)'^ can be compufed in 
time if two m x m Boolean matrices can be multiplied in time 0{nrG). 

Proof 

The result of fhe algorithm above is guaranteed to be the same as the result of matrix 
multiplication using the 0 operation of Figure 5 because if considers all combinations 
of z, j, and k and all pairs of nonferminals and copy symbols, and applies fhe same sef 
of conditions. This is possible because each of the conditions in Figure 5 applies either 
to a pair (z, k) or (fc, j), in which case we apply the condition to input matrices to the 
Boolean matrix multiplication, or to the pair (z, j), in which case we apply the condition 
to the result of the Boolean matrix multiplication. Crucially, no condition in Figure 5 
involves z, j, and k simultaneously. 

The Boolean matrix algorithm takes time for each mafrix multiplication, 

while the pre- and post-processing steps for each matrix multiplication take only 
The number of Boolean mafrix mulfiplicafions depends on the grammar, but 
is constant with respect to rz, yielding an overall runtime of OirG'^). ■ 

The final parsing algorifhm is given in Figure 8. If works by compufing the seed 
matrix T, and then finding its transitive closure. Finally, it checks whether the start 
symbol appears in a cell with an address that spans the whole string. If so, fhe sfring is 
in the language of the grammar. 
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Inputs: An LCFRS grammar as defined in §2 that satisfies Condition 4.1 and a sentence 

Wl • • • Wn- 

Outputs: True if wi • • • Wn is in the language of the grammar. False otherwise. 

Algorithm: 

1. Compute T as the seed matrix using the algorithm in Figure 4. 

2. Compute the transitive closure of T with the multiplication operator in Figure 5 
and using Boolean matrix multiplication (§4.5). 

3. Return True if {S, (0), (n)) belongs to the cell ((0), (n)) in the computed transitive 
closure, and False otherwise. 


Figure 8 

Algorithm for recognizing binary linear context-free rewriting systems when Condition 4.1 is 
satisfied by fhe LCFRS. 


5. Computational Complexity Analysis 

As mentioned in the previous section, the algorithm in Figure 8 finds the transitive 
closure of a matrix under our definition of matrix multiplication. The operations U and 
C) used in our matrix multiplication distribute. The C) operator takes the cross product 
of two sets, and applies a filtering condition to the results; the fact that (x C) j/) U (x C) 
z) = X ^ {y U x) follows from the fact that it does not matter whether we take the cross 
product of the union, or the union of the cross product. However, unlike in the case of 
standard matrix multiplication, our (g) operation is not associative. In general, x (g (y (g 
2 :) yf (x (g y) g) z, because the combination of y and z may be allowed by the LCFRS 
grammar, while the combination of x and y is not. 

Lemma 8 

The transitive closure of a matrix of size (2n)‘^ x (2n)'^ can be computed in time 
if 2 < w < 3, and two m x m Boolean matrices can be multiplied in time 0{m‘^). 

Proof 

We can use the algorithm of Valiant for finding the closure of upper triangular matrices 
under distributive, non-associative matrix multiplication. Because we can perform one 
matrix product in time by Lemma 7, the algorithm of Valiant (1975, Theorem 2) 

can be used to compute transitive closure also in time 0 {n‘^‘^). ■ 

When Valiant's paper was published, the best well-known algorithm known for 
such multiplication was Strassen's algorithm, with M{n) = Since then, it is 

known that M{n) = 0(n‘^) for uj < 2.38 (see also §1). There are ongoing attempts to 
further reduce w, or find lower bounds for M(n). 

The algorithm for transitive closure gives one of the main results of this article: 

Theorem 1 

A single-initial binary LCFRS meeting Condition 4.1 can be parsed in time 
where d is the contact rank of the grammar, 2 < w < 3, and two m x m Boolean matrices 
can be multiplied in time 0 {m‘^). 
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Proof 

By Lemma 8, step 2 of the algorithm in Figure 8 takes By Lemma 5, the result of 

step 2 gives all nonterminals that are derivable for the given spans ■ 

Parsing a binary LCFRS rule with standard chart parsing techniques requires time 
0 (^n'^(^)+v(B)+<p(C)y p _ ceK + ‘fiiC)). The worst-case com¬ 

plexity of LCFRS chart parsing techniques is 0(riF). We can now ask the question: in 
which case the algorithm in Figure 8 is asymptotically more efficient than standard chart 
parsing techniques with respect to n? That is, in which cases is = o^nP)! 

Clearly, this would hold whenever duj < p. By definition of d and p, a sufficient 
condition for that is that for any rule BC & TZ it holds that:^ 

f ‘P{A) + - Lp{C), 'I , 

max<^ p{A) - ip{B) + p{C), \ < - {p{A) + p{B) + p(C)). 
{-p{A)+p{B)+p{C)] ^ 

This means that for any rule, the following conditions should hold: 

uj{if{A) + (f{B) - p(C)) < <p{A) + (p{B) + 

~ viB) + p{C)) < p{A) + (p{B) + ip{C), 

Lo{—ip{A) -\- (p{B) + (p{C)) < (p{A) + (p{B) + (p{C). 

Algebraic manipulation shows that this is equivalent to having: 

(p{A) + <p(B) < ‘f(C'), 

p{B) + (p{C) < 

(p(C) + ip{A) < ‘f(-S)- 


For the best well-known algorithm for matrix multiplication, it holds that: 

w -I- 1 


w — 1 


> 2.44. 


For Strassen's algorithm, it holds that: 

w -I- 1 

(jj — 1 


> 2.06. 


We turn now to analyze the complexity of the algorithm in Figure 9, giving the main 
result of this article for arbitrary LCFRS: 


4 For two sets of real numbers, X and Y, it holds that if for all a; G X there is a j/ G F such that x < y, then 
maxX < maxY. 
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Inputs: An LCFRS grammar as defined in §2 and a sentence wi ■ ■ -Wn- 
Outputs: True if wi • • • Wn is in the language of the grammar. False otherwise. 

Algorithm: 

1. Compute T as the seed matrix using the algorithm in Figure 4. 

2. Repeat until T does not change: T -i— (n(T))^. 

3. Return True if {S, (0), (n)) belongs to the cell ((0), (n)) in the computed transitive 
closure, and False otherwise. 


Figure 9 

Algorithm for recognizing binary LCFRS when Condition 4.1 is not necessarily satisfied by the 
LCFRS. 


Theorem 2 

A single-initial binary LCFRS can be parsed in time where d is the contact 

rank of the grammar, 2 < w < 3, and two m x m Boolean matrices can be multiplied in 
time 0 (to“). 


Proof 

The algorithm of Figure 9 works by iteratively applying the transitive closure and the 
copying operator until convergence. At convergence, we have recognized all derivable 
nonterminals by Lemma 6. Each transitive closure has the asymptotic complexity of 
by Lemma 8. Each It application has the asymptotic complexity of 0{n‘^). As 
such, the total complexity is where t is the number of iterations required to 

converge. At each iteration, we discover at least one new nonterminal. The total number 
of nodes in the derivation for the recognized string is 0(n) (assuming no unary cycles 
or e rules). As such t = 0(n), and the total complexity of this algorithm is ■ 


6. Applications 

Our algorithm is a recognition algorithm which is applicable to binary LCERS. As such, 
our algorithm can be applied to any LCERS, by first reducing it to a binary LCERS. We 
discuss results for specific classes of LCERS in this section, and return to the general 
binarization process in §7.6. 

LCERS subsumes context-free grammars, which was the formalism that Valiant 
(1975) focused on. Valiant showed that the problem of CEG recognition can be reduced 
to the problem of matrix multiplication, and as such, the complexity of CEG recognition 
in that case is 0(n“). Our result generalizes Valiant's result. CEGs (in Chomsky normal 
form) can be reduced to a binary LCERS with / = 1. As such, d = 1 for CEGs, and our 
algorithm yields a complexity of 0(n“). (Note that CEGs satisfy Condition 4.1, and 
therefore we can use a single transitive closure step.) 

LCERS is a broad family of grammars, and it subsumes many other well-known 
grammar formalisms, some of which were discovered or developed independently of 
LCERS. Two such formalisms are tree-adjoining grammars (Joshi and Schabes 1997) 
and synchronous context-free grammars. In the next two sections, we explain how our 
algorithmic result applies to these two formalisms. 
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6.1 Mildly Context-Sensitive Language Recognition 

Linear context-free rewriting systems fall under the realm of mildly context-sensitive 
grammar formalisms. They subsume four important mildly context-sensitive for¬ 
malisms that were developed independently and later shown to be weakly equivalent 
by Vijay-Shanker and Weir (1994): tree-adjoining grammars (Joshi and Schabes 1997), 
linear indexed grammars (Gazdar 1988), head grammars (Pollard 1984) and combina¬ 
tory categorial grammars (Steedman 2000). Weak equivalence here refers to the idea 
that any language generated by a grammar in one of these formalisms can be also be 
generated by some grammar in any of the other formalisms among the four. It can be 
verified that all of these formalisms are unbalanced, single-initial LCFRSs, and as such, 
the algorithm in Figure 8 applies to them. 

Rajasekaran and Yooseph (1998) showed that tree-adjoining grammars can be 
parsed with an asymptotic complexity of 0{M(n^)) = 0(n^'^®). While he did not dis¬ 
cuss that, the weak equivalence between the four formalisms mentioned above implies 
that all of them can be parsed in time 0{M{n?)). Our algorithm generalizes this result. 
We now give the details. 

Our starting point for this discussion is head grammars. Head grammars are a 
specific case of linear context-free rewriting systems, not just in the formal languages 
they define - but also in the way these grammars are described. They are described 
using concatenation production rules and wrapping production rules, which are di¬ 
rectly transferable to LCFRS notation. Their fan-out is 2. We focus in this discussion 
on “binary head grammars," defined analogously to binary LCFRS - the rank of all 
production rules has to be 2. The contact rank of binary head grammars is 2. As such, 
our paper shows that the complexity of recognizing binary head grammar languages is 
0(M(n^)) = 0(n‘*’^®). 

Vijay-Shanker and Weir (1994) show that linear indexed grammars (LIGs) can ac¬ 
tually be reduced to binary head grammars. Linear indexed grammars are extensions 
of CFGs, a linguistically-motivated restricted version of indexed grammars, the latter 
of which were developed by Aho (1968) for the goal of handling variable binding 
in programming languages. The main difference between LIGs and CFGs is that the 
nonterminals carry a “stack," with a separate set of stack symbols. Production rules 
with LIGs copy the stack on the left-hand side to one of the nonterminal stacks in the 
righthand side,® potentially pushing or popping one symbol in the new copy of the 
stack. For our discussion, the main important detail about the reduction of LIGs to head 
grammars is that it preserves the rank of the production rules. As such, our paper shows 
that binary LIGs can also be recognized in time 0(n‘*’^®). 

Vijay-Shanker and Weir (1994) additionally address the issue of reducing combina¬ 
tory categorial grammars to LIGs. The combinators they allow are function application 
and function composition. The key detail here is that their reduction of CCG is to an 
LIG with rank 2, and as such, our algorithm applies to CCGs as well, which can be 
recognized in time 0(n'^'^®). 

Finally, Vijay-Shanker and Weir (1994) reduced tree-adjoining grammars to combi¬ 
natory categorial grammars. The TAGs they tackle are in “normal form," such that the 
auxiliary trees are binary (all TAGs can be reduced to normal form TAGs). Such TAGs 
can be converted to weakly equivalent GCG (but not necessarily strongly equivalent), 
and as such, our algorithm applies to TAGs as well. As mentioned above, this finding 


5 General indexed grammars copy the stack to multiple nonterminals on the right-hand side. 
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supports the finding of Rajasekaran and Yooseph (1998), who showed that TAG can be 
recognized in time 0(M(n^)). 

For an earlier discussion connections between TAG parsing and Boolean matrix 
multiplication, see Satta (1994). 

6.2 Synchronous Context-Free Grammars 

S 5 mchronous Context-Free Grammars (SCFGs) are widely used in machine translation 
to model the simultaneous derivation of translationally equivalent strings in two natu¬ 
ral languages, and are equivalent to the S 5 mtax-Directed Translation Schemata of Aho 
and Ullman (1969). SCFGs are a subclass of LCFRS where each nonterminal has fan¬ 
out two: one span in one language and one span in the other. Because the first span of 
the l.h.s. nonterminal always contains spans from both r.h.s. nonterminals, SCFGs are 
always single-initial. Binary SCFGs, also known as Inversion Transduction Grammars 
(ITGs), have no more than two nonterminals on the r.h.s. of a rule, and are the most 
widely used model in syntax-based statistical machine translation. 

S 5 mchronous parsing with traditional tabular methods for ITG is 0(n®), as each of 
the three nonterminals in a rule has fan-out of two. ITGs, unfortunately, do not satisfy 
Condition 4.1, and therefore we have to use the algorithm in Figure 9. Still, just like 
with TAG, each rule combines two nonterminals of fan-out two using two combination 
points. Thus, d = 2, and we achieve a bound of for ITG, which is 

using the current state of the art for matrix multiplication. 

We achieve even greater gains for the case of multi-language s 5 mchronous parsing. 
Generalizing ITG to allow two nonterminals on the righthand side of a rule in each 
of k languages, we have an LCFRS with fan-out k. Traditional tabular parsing has an 
asymptotic complexity of 0{n^^), while our algorithm has the complexity of ©(n^^"*"^). 

Another interesting case of a synchronous formalism that our algorithm improves 
the best-well known result for is that of binary synchronous TAGs (Shieber and Schabes 
1990) - i.e. a TAG in which all auxiliary trees are binary. This formalism can be reduced 
to a binary LCFRS. A tabular algorithm for such grammar has the asymptotic complex¬ 
ity of 0{n}^). With our algorithm, d = 4 for this formalism, and as such its asymptotic 
complexity in that case is 0(n® ®^). 

7. Discussion and Open Problems 

In this section, we discuss some extensions to our algorithm and open problems. 

7.1 Turning Recognition into Parsing 

The algorithm we presented focuses on recognition: given a string and a grammar, 
it can decide whether the string is in the language of the grammar or not. From an 
application perspective, perhaps a more interesting algorithm is one that returns an 
actual derivation tree, if it identifies that the string is in the language. 

It is not difficult to adapt our algorithm to return such a parse, without changing 
the asymptotic complexity of Once the transitive closure of T is computed, 

we can backtrack to find such parse, starting with the start symbol in a cell spanning the 
whole string. When we are in a specific cell, we check all possible combination points 
(there are d of those) and nonterminals, and if we find such pairs of combination points 
and nonterminals that are valid in the chart, then we backtrack to the corresponding 
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Figure 10 

Upper left: Combination of spans for SCFG rule [S ABC D, B D AC]. Upper right and 
bottom row: three steps in parsing binarized rule. 


cells. The asymptotic complexity of this post-processing step is which is less 

than (uj > 2, d > 1). 

This post-processing step corresponds to an algorithm that finds a parse tree, given a 
pre-calculated chart. If the chart was not already available when our algorithm finishes, 
fhe asymptotic complexity of this step would correspond to the asymptotic complexity 
of a naive tabular parsing algorithm. It remains an open problem to adapt our algo¬ 
rithm to probabilistic parsing, for example - finding the highest scoring parse given a 
probabilistic or a weighted LCFRS (Kallmeyer and Maier 2010). See more details in §7.3. 

7.2 General Recognition for Synchronous Parsing 

Similarly to LCFRS, the rank of an SCFG is the maximal number of nonterminals that 
appear in the right-hand side of a rule. Any SCFG can be binarized into an LCFRS 
grammar. However, when the SCFG rank is arbitrary, this means that the fan-out of 
fhe LCFRS grammar can be larger than 2. This happens because binarization creates 
intermediate nonterminals that span several substrings, denoting binarization steps of 
fhe rule. These substrings are eventually combined into two spans, to yield the language 
of the SCFG grammar (Huang et al. 2009). 

Our algorithm does not always improve the asymptotic complexity of SCFG pars¬ 
ing over tabular methods. For example. Figure 10 shows the combination of spans for 
the rule [S' AB CD,B D AC], along with a binarization into three simpler LCFRS 
rules. A naive tabular algorithm for this rule would have the asymptotic complexity of 
0(n^°), but the binarization shown in Figure 10 reduces this to 0(n®). Our algorithm 
gives a complexity of 0(n®-®^), as the second step in the binarization shown consists of 
a rule with d = 4. 

7.3 Generalization to Weighted Logic Programs 

Weighted logic programs (WLPs) are declarative programs, in the form of Horn clauses 
similar to those that Prolog uses, that can be used to formulate parsing algorithms 
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such as CKY and other types of dynamic programming algorithms or NLP inference 
algorithms (Eisner, Goldlust, and Smith 2005; Cohen, Simmons, and Smith 2011). 

For a given Horn clause, WLPs also require a "join" operation that sums (in some 
semiring) over a set of possible values in fhe free variables in the Horn clauses. With 
CKY, for example, this sum will be performed on fhe mid-point concatenating two 
spans. This join operation is also the t 5 rpe of operation we address in this paper (for 
LCFRS) in order to improve their asymptotic complexity. 

It remains an open question to see whether we can generalize our algorithm to 
arbitrary weighted logic programs. In order to create an algorithm that takes as input 
a weighted logic program (and a set of axioms) and "recognizes" whefher the goal is 
achievable, we would need to have a generic way of specifying the set N, which was 
specialized to LCFRS in this case. Not only that, we would have to specify N in such a 
way fhat fhe asymptotic complexity of fhe WLP would improve over a simple dynamic 
programming algorifhm (or a memoizafion technique). 

In addition, in this paper we focus on the problem of recognition and parsing for 
unweighted grammars. Benedi and Sanchez (2007) showed how to generalize Valiant's 
algorithm in order to compute inside probabilities for a PCFG and a string. Even if we 
were able to generalize our addressing scheme to WLPs, it remains an open question to 
see whether we can go beyond recognition (or unweighted parsing). 

7.4 Rytter's Algorithm 

Rytter (1995) gives an algorithm for CFG parsing with the same time complexity 
as Valiant's, but a somewhat simpler divide-and-conquer strategy. Rytter's algorithm 
works by first recursively finding all charf items entirely within the first half of fhe string 
and entirely within the second half of the string. The combination step uses a shortest 
path computation to identify the sequence of chart items along a spine of the final parse 
tree, where the spine extends from fhe root of fhe tree to the terminal in position n/2. 
Rytter's algorithm relies on the fact that this spine, consisting of chart items that cross 
the midpoint of fhe string, forms a single path from the root to one leaf of the derivation 
tree. This property does not hold for general LCFRS, since two siblings in the derivation 
tree may both correspond to multiple spans in the string, each containing material on 
both sides of the string midpoint. For this reason, Rytter's algorithm does not appear to 
generalize easily to LCFRS. 

7.5 Relation to Multiple Context-Free Grammars 

Nakanishi et al. (1998) developed a matrix multiplication parsing algorithm for multiple 
context-free grammars (MCFGs). When these grammars are given in a binary form, 
they can be reduced to binary LCFRS. Similarly, binary LCFRS can be reduced to binary 
MCFGs. The algorithm that Nakanishi et al. develop is simpler than ours, and does not 
directly tackle the problem of fransitive closure for LCFRS. More specifically, Nakanishi 
et al. multiply a seed matrix such as our T by itself in several steps, and then follow up 
with a copying operation between equivalent cells. They repeat this n times, where n is 
the sentence length. As such, the asymptotic complexity of fheir algorithm is identical 
for both balanced and unbalanced grammars, a distinction they do not make. 

The complexity analysis of Nakanishi et al. is different than ours, but in certain 
cases, yields identical results. For example, if (p(a) = / for all a € C, and the grammar is 
balanced, then both our algorithm and their algorithm give a complexity of 
If fhe grammar is unbalanced, fhen our algorithm gives a complexity of 0{n‘^^), while 
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the asymptotic complexity of their algorithm remains As such, Nakanishi et 

al.'s algorithm does not generalize Valiant's algorithm - its asymptotic complexity for 
context-free grammars is and not 0{n‘^). 

Nakanishi et al. pose in their paper an open problem, which loosely can be re¬ 
worded as the problem of finding an algorithm that computes the transitive closure 
of T without the extra Oin) factor that their algorithm incurs. In our paper, we provide 
a solution to this open problem for the case of single-initial, unbalanced grammars. The 
core of the solution lies in the matrix multiplication copying mechanism described in 
§4.2. 

7.6 Optimal Binarization Strategies 

The two main grammar parameters that affect the asymptotic complexity of parsing 
with LCFRS (in its general form) are the fan-out of the nonterminals and the rank of 
the rules. With tabular parsing, we can actually refer to the parsing complexity of a 
specific rule in the grammar. Its complexity is 0{vP), where the parsing complexity p is 
the total fan-out of all nonterminals in the rule. For binary rules of the form A ^ BC, 
P = + V{B) +^{0). 

To optimize the tabular algorithm time complexity of parsing with a binary LCFRS, 
equivalent to another non-binary LCFRS, we would want to minimize the time com¬ 
plexity it takes to parse each rule. As such, our goal is to minimize p^A) + p{B) + (p{C) 
in the resulting binary grammar. Gildea (2011) has shown that this metric corresponds 
to the tree width of a dependency graph which is constructed from the grammar. It is 
not known whether finding the optimal binarization of an LCFRS is an NP-complete 
problem, but Gildea (2011) shows that a polynomial time algorithm would imply im¬ 
proved approximation algorithms for the treewidth of general graphs. 

In general, the optimal binarization for tabular parsing may not by the same as the 
optimal binarization for parsing with our algorithm based on matrix multiplication. In 
order to optimize the complexity of our algorithm, we want to minimize d, which is the 
maximum over all rules A ^ B C of 

d{A BC) = max{(/?(A) -|- p{B) — cp{C), <p(A) — cp{B) + (f{C), —<p(A) + <p{B) + 

For a fixed binarized grammar, d is always less than p, the tabular parsing com¬ 
plexity, and, hence, the optimal d* over bmarizations of an LCFRS is always less than 
the optimal p* for tabular parsing. However, whether any savings can be achieved with 
our algorithm depends on whether cod* < p*, or cod* -I-1 < p* in the case of balanced 
grammars. Our criterion does not seem to correspond closely to a well-studied graph- 
theoretic concept such a treewidth, and it remains an open problem to find an efficient 
algorithm that minimizes this definition of parsing complexity. 

It is worth noting that d{A BC) + p{B) + p(G)). As such, this gives a 

lower bound on the time complexity of our algorithm relative to tabular parsing using 
the same binarized grammar. If 0{n}^ ) is the asymptotic complexity of our algorithm, 

and 0{n*^) is the asymptotic complexity of a tabular algorithm, then — > — > 0.79. 

G 3 


8. Conclusion 

We described a parsing algorithm for binary linear context-free rewriting systems that 
has the asymptotic complexity of 0(n‘^'^+^) where oj < 2.38, d is the "contact rank" of 
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Symbol 

Description 

1st mention 

M(n) 

The complexity of Boolean n x n matrix multiplication 

§1 

CO 

Best well-known complexity for M{n), M{n) = 0{'nA) 

§1 

[n] 

Set of integers {1,..., n} 

§2 

[n]o 

[n] U {0} 

§2 

C 

Nonterminals of the LCPRS 

§2 

r 

Terminal symbols of the LCPRS 

§2 

V 

Variables that denote spans in grammar 

§2 

n 

Rules in the LCFRS 

§2 

A,B,C 

Nonterminals 

§2 

f 

Maximal fan-out of the LCFRS 

Eq.3 

ip{A) 

Fan-out of nonterminal A 

§2 

y 

Denoting a variable in V (potentially subscripted) 

§2 

T 

Seed matrix 

§3 

N, N{d) 

Set of indices for addresses in the matrix 

Eq.4 

Cj 

Indices for cells in T. i, j G N 

§4.1 

d 

Grammar contact rank 

§4.1 

M 

Tij is a subset of M 

§4.1 


Copying/marking symbols for rows 

§4.1 


Copying/marking symbols for columns 

§4.1 

n 

Length of sentence to be parsed 

§1 

< 

Total order between the set of indices of T 

§4.1 


Merged sorted sequence of i U j, divided into pairs 

§4.1 

remove (t!, x) 

Removal of x from a sequence v 

Figure 5 

insert(v, x) 

Insertion of a; in a sequence v 

§4.5 

n 

Copying operator 

§4.5 


Table 1 

Table of notation symbols used in this paper. 


the grammar (the maximal number of combination points in the rules in the grammar 
in single-initial form) and n is the string length. Our algorithm has the asymptotic 
complexity of for a subset of binary LCFRS which are unbalanced. Our result 

generalizes the algorithm of Valiant (1975), and also reinforces existing results about 
mildly context-sensitive parsing for tree-adjoining grammars (Rajasekaran and Yooseph 
1998). Our result also implies that inversion transduction grammars can be parsed 
in time and that synchronous parsing with k languages has the asymptotic 

complexity of where k is the number of languages. 

A. Notation 

Table 1 gives a table of notation for symbols used throughout this paper. 
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